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I. INTRODUCTION 



In 2003 BaBar Collaboration discovered a positive-parity scalar charm strange meson 
Dsj(2317) with a very narrow width which was confirmed by CLEO j3| later. In the 
same experiment CLEO also observed the 1+ partner state at 2460 MeV jl]. Since these two 
states lie below the DK and D*K threshold, respectively, the potentially dominant s-wave 
decay modes Dsj(2317) DK etc., are kinematically forbidden. Thus the radiative decays 
and isospin-violating strong decays become the dominant decay modes. Therefore both of 
them are very narrow. 

The discovery of these two states has triggered heated discussion on their nature in 
literature. The key point is to understand their low masses. The mass of i5sj(2317) is 
significantly lower than the expected values in the range of 2.4 — 2.6 GeV within quark 
models j^. The model using the heavy-quark mass expansion of the relativistic Bethe- 
Salpeter equation predicted a lower value 2.369 GeV which is still higher than the 
experimental data by about 50 MeV. 

From the experience with ao//o(980). Van Beveren and Rupp argued that the low mass 
of Dsj(2317) could arise from the mixing between the 0+ cs state and the DK continuum. 
In this way the lowest O"*" state could be pushed much lower than that expected from the 
quark models. 

The mass of Ds{Q^) state from the lattice QCD calculation is also significantly larger than 
the experimentally observed mass of Dsj(2317) [sl, 0, S]- It is also pointed out in Ref. 
that Dsj(2317) might receive a large component of DK continuum, which makes the lattice 
simulation very difficult. 

The difficulty with the cs interpretation leads many authors to speculate that Dsj{2ZVI) 
is a cqsq four quark state [ol, [l^, or a strong D-n atom [11]. However, calculations based 
on the quark model show that the mass of the four quark state is much larger than that 
of the 0+ cs state [ll, [iHl- The radiative decay of Z)sj(2317) also favors that it is a cs 
state IJ]. Furthermore, there are two 0'^ states in the four quark system and one in the 
two-quark system. Only one 0^ state has been found below the 2.86 GeV resonance in the 
experimental search by BaBar [15], consistent with the cs interpretation. 

This pro blem has been treated with QCD sum rules in the heavy quark effective theory 
in Ref. [l6|. The resulting Ds(O^) mass is consistent with the experimental data within 
large theoretical uncertainties. However, the central value is still larger than the data by 
90 MeV. Even larger result for the Ds(0"'") mass was obtained in the earlier work with the 
sum rule in full QCD V]\. It has been pointed out in Ref. 16| that, in the formalism of 
QCD sum rules, the physics of mixing with DK continuum resides in the contribution of 
DK continuum in the sum rule, and including this characteristic contribution should render 
the mass of -Ds(O^) lower. 

Recently, there have been two investigations on this problem using sum rules in full 
QCD including the 0{as) corrections. In Ref. 18| the value of the charm quark pole mass 
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Mc = 1.46 GeV is used, and the mass of 0^ cs state is found to be 100 — 200 MeV higher 
than the experimental data. On the other hand, in Ref. [1^] the current quark mass rric = 
1.15 GeV (corresponding to Mc ~ 1.3 GeV to 0{as)) is used, and the central value of the 
resulting O"*" cs mass is in agreement with the data. However, a low value of rric is used and 
the same value of the continuum threshold (denoted by sq below) is used for O"*" cs and 0~ 
cs. 

On the other hand, the perturbative three loop, order correction to the two-point 



correlation function with one heavy and one massless quark has been calculated 20|, l2l| . It 



turns out that in the pole mass scheme used by many previous analyses including Ref. [l8|, 



191], the perturbative expansion is far from converging. However, taking the quark mass in 



the modified minimal subtraction (MS) scheme 22|], better convergence of the higher order 



corrections is obtained, and thus a more reliable determination of physical quantities of the 



lowest lying resonances becomes feasible 23 ] 



Usually the contribution of two-particle continuum is neglected within the QCD sum rule 
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25 and the 



formalism. However, because of the large s-wave coupling of Ds{0~^)DK 
adjacency of the Ds{0^) mass to the DK threshold, this contribution may not be neglected 
in the considered case. In the present work, we shall therefore calculate this contribution 
and include it in the QCD sum rule. In the meantime we take into account the perturbative 
three loop order correction and work in the MS scheme. We find that the DK continuum 
contribution indeed renders both the mass and the decay constant of 1)^(0"'") significantly 
lower. 

In Section [ITl we give a short overview of the traditional QCD sum rule for the scalar 
charm strange meson. Then we derive the DK continuum contribution and write down the 
full sum rule in Section IIIII Finally, the numerical results and our conclusions are presented 
in Section [IVl Some relevant formulas and expressions used in this paper are collected in 
Appendices. 

II. THE TRADITIONAL QCD SUM RULE FOR THE SCALAR CHARM- 
STRANGE MESON 

We consider the scalar correlation function 

U{p') ^ z ldxe'P^O\T{j{x)j{0)^\ 0), (1) 

where the renormalization invariant operator j{x) is defined as 

j(x) = (mc - nis) : s(x) c(x) : , (2) 

with rric and being the charm and strange quark current mass, respectively. Up to a 
subtraction polynomial in j9^, the correlation function H(p^) satisfies the following dispersion 
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relation 



Il(p^) = I ds + subtractions . (3) 

' [s — — le) 







At the quark gluon level, the spectra function p{s) is calculable using the renormaliza- 
tion group improved perturbation theory in the framework of the operator product expan- 



sion (OPE). Following Jamin and Lange 23|], in this paper we shall adopt the MS running 
quark mass scheme rather than the pole mass one, and take into account both the 0{a^) 
terms in the perturbation theory obtained in Ref. 0, El| and the corrections from the 
light quark mass up to order m^. In addition we have included the contribution from the 
four quark condensation which affects the final result for Ds{0'^) mass only by a few Mev. 
For convenience, all the relevant expressions for p{s) at the quark-gluon level, denoted by 
Pqct){s), are summarized in Appendix A. 

On the other hand, p{s) can be phenomenologically written in terms of contributions from 
intermediate hadronic states. Generally, the spectral density at the hadronic level, denoted 
by ph, is taken to be the pole term of the lowest lying hadronic state plus the continuum 
starting from some threshold, with the latter identified with the QCD continuum 

^"^^^ - foM^6{t - M^) + QCD continuum x e{t - sq) , (4) 



TT 



where /o is the vector current decay constant of 0"*" cs particle, analogous to = 131 MeV. 
M is the mass of this particle, and Sq is the continuum threshold above which the hadronic 



spectral density is modeled by that at the quark gluon level. The recent works [18|, Il9[ also 
use the above ansatz. 

After making the Borel transformation to suppress the contribution of higher excited 
states and invoking the quark-hadron duality, one arrives at the sum rule 

oo 

dt exDf-— 1 - fdt ^Q^"^^^ exDf-— 1 (5) 

dt ^ exp[ ^2J-y ^ ^2 J- (5) 

M2 



Following Ref. [23j, the lower limit of the integration in the above equation is taken to be the 
charm quark pole mass Mc, which can be expressed in terms of the running mass mc{pm) 
through the perturbative three-loop relation as defined in Appendix B. 



III. THE CONTRIBUTION OF DK CONTINUUM 

The contribution of two-particle continuum to the spectral density can safely be neglected 
in many cases, as usually done in the traditional QCD sum rule analysis. One typical example 
is the p meson sum rule, where the two pion continuum is of p-wave nature. Its contribution 
to the spectral density is tiny and the p pole contribution dominates. 
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However, there may be an exception when the O"*" particle couples strongly to the two- 
particle continuum via s-wave. In such case, there is no threshold suppression and the 
two-particle continuum contribution may be more significant. The strong coupling of the O"*" 
particle with the two-particle state and the adjacency of the O"*" mass to the DK continuum 
threshold result in large coupling channel effect, which corresponds to the configuration of 
mixing in the formalism of quark model. In the problem under consideration, the mass 
of -Dsj (231 7) is only about 45 MeV below the DK threshold, and the s-wave coupling of 
Da{Q^)DK is found to be very large 2J, |25|]. Therefore, one may have to take into account 
the DK continuum contribution carefully. 

The importance of Dvr continuum contribution in the sum rule for -D(O^) meson was first 
emphasized in Ref. [2^, based on the duality consideration in the case where the -D(O^) 
mass is higher than the Dt[ threshold. Based on the soft pion theorem, two of us also made 
a crude analysis of the i?7r continuum contribution in the case where the O"*" particle mass is 



higher than the threshold [25j. In this work, we calculate the continuum contribution more 
carefully in the case where the O"*" particle mass is lower and very close to the two-particle 
continuum threshold. 

Let F{t) be the form-factor defined by 

F(t) = (0|c(0).(0)|M). (6) 

From the large s-wave coupling of Ds{{)^)DK and the adjacency of the Ds(0"'") mass to the 
DK threshold, one expects that in the low energy region, F{t) is dominated by the product 
of a factor of the Ds{{)^) pole and a factor from the final state interactions. In the low energy 
region with {niD + itikY < t < sq < 8 GeV^ needed in our sum rule, the effect of inelastic 
DK scattering is suppressed by the phase space. Therefore, we take the approximation to 
consider only the DK scattering with only elastic intermediate states. It can be described 
by the Ds{0~^)DK interaction and the DDKK chiral interaction in the low energy effective 
lagrangian, which can be represented by a series of bubble diagrams shown in Fig. [H 



o 





FIG. 1: Heavy, light, and dotted lines represent Z)s(O^), D, and K, respectively. Black circle 
represents the Born s-wave amplitude of DK scattering, and blank one the scalar current. 



The s-wave Born amplitude of DK scattering represented by the black circles in Fig. [T] 
contains three terms. The first one is the t channel pole term jzm^i with qq being the 
Ds{Q^)DK coupling constant and Mq being the mass parameter normalized at the scale m|) 
in the effective lagrangian. 
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The second term corresponds to the direct DDKK interaction in the effective lagrangian. 
Let p, k and p', k' be the four momentum of D, K mesons in the initial and final state 
respectively, and s = {p — k'Y = {p' — k)"^. In the chiral effective lagrangian, the amplitudes 
for the processes D+K'^ D+K°, D^K+ D°K+, and D+K'^ ^ D^K+ in the low energy 
ko region of K meson needed in the QCD sum rule are all equal to 



+1 

-1 



2Vt{kQ + fco) - 2kok'Q -k 



12 



(7) 



in the center of mass system. Here k and k' should be separately included in the integrals 
for two adjacent loops. 

The third term of the s-wave Born amplitude arises from the crossing s channel pole term 



■9l 



+1 



dcosO 



B 



.91 1 , A 

I— — In — 

2 B A + B 



(8) 



where 



A = t-2Vt{ko + k'Q) +2kok'Q + 2mj^- M^, B = 2\k\\k'\ . (9) 

For simplicity, we put the on-shell values of ko, /cq, \k\, \k'\ into A and B in the above 
equations. The effect of this approximation on our final results is expected to be small, 
since the contribution of the s channel pole term is relatively small. As a result, one finds 
that the s channel pole term is an analytic function of t with only a short cut, the length of 
which is only 0.146 GeV^ for experimental values of the corresponding masses. Therefore, 

2 

it can be well approximated by a pole form —i-^zf^, where 



1 
2 



2m1 + 2m^ 



K 



[m 



D 



m 



2 

k) 



(10) 



2m% + 2ml - - 



K 



to 



2ml 



2ml 



(11) 



With the above results for the three terms of the s-wave Born amplitude, we can now 
evaluate the sum of the series of the bubble diagrams shown in Fig. [TJ Let fn{t) be the 
partial sum of the series of loop diagrams in Fig. [H with the loop number less or equal to 
n. It can then be written in the form 



-1 



2/ 



\2Vik',-k'')-2l 



9l 



t 



Mi 



+ 



t-i 







fnO + 2{Vt-k',)U-fn2\, (12) 



where k^ and k' is the energy and momentum of final-state K meson, and the three unknown 
functions /„« {i = 0, 1, 2) correspond to the diagrams with a factor 1, ko, and k"^ respectively 
at the last vertex, which contributes to the integration over the last loop of each diagram. 
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Let be integrals defined by Eqs.(Cl)-(C6) which appear as the loop integrals of 

the individual loop diagrams shown in Fig. [TJ They can be evaluated using dimensional 
regularization 27|, with the corresponding analytic results given in (C1)-(C6). A recurrence 
relations can be written between /„(t) and and hence between fniit) and /(n+i)i(t), 

the coefficients of which are linear combinations of the loop-integral functions Taking 
the limit lim„^oo/m(^) = fi(t), lim„^oo/n(^) = fit) and separating out terms with the 
factor 1, /cq, and k''^, we can obtain a system of three linear equations for the three unknown 
functions fi(t) 

h (2E2 + E3) + 2/oS4] + 4 (/1E1 + /0E2) t 

c 1 



2 /i E4 + /o E5 — 2 

+4 gl L So - /o [Ml + gl Eo - t) 



it - tof (t - M2) (t - to) 



go (l-2/o^7o) (2EiVt-E3) +2/1 (^o' Si - (M^ + (?o' So) Vi + n 



t - A'Q 



+ /2 



+2f 



E3-2EiV^+2/i,,^Eo , 
2cgl [/iEi + /oS3-(/iSo + 2/oSi) 



2/ 



+4c/o/^(7o'S 



t-tn 



+ 



(t-to)2 (t-M2) (t-to) 







(13) 



/2S1 + /0S4- [/2S0 + /0 (2S2 + S3)] Vt + 2/i (-/|: + S2-2Siv^+Sot 

00 (1 - /o ^0) Si + k M2 - ^0 (1 - /o ^0) Sol Vt-fot'^ 
+'^fK 7rr— + 2/oSit 



+2/ 



cfogl (Si-Sov^) 



K 



t-to 

2/1S1 + /0 



0, 



(14) 



S3-2SiVt + 2/1-^72 So f 



2/^ 



(15) 



from which the analytic forms for fi{t) can then be deduced. The results are shown in Eqs. 
(C9)-(C12). 

Finally, with the explicit expressions for fi{t) and Sj(t) given in Appendix C, and putting 
the on-shell value of k' to Eq. (|T2l) . we obtain 



-9o 



t - Ml 2fl 
A 

t-M^- A{t) 



c9o 



t-M^ t-to 



fo + 



/1-/2 



(16) 
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where A(t) is given by Eq. (IClSp . 

Similarly, with the same series of DK loops included, the full propagator of Ds(0"*") meson 
is related to the function /o(t) through 

Prop(t) = [1 - gMt)] . (17) 

Using the solution for /o(t) obtained above and given by Eq. (]C9I) . it can be further rewritten 



as 

with Ai(t) given by Eq. flUTil) . 

We have chosen the scale /i so that the mass parameter Mq is the physical mass of Ds{0^) 
in our approximation, i.e., Ai(Mq) = 0. The bare coupling constant is related to the 
physical coupling constant g hy g = g^/y/Z, where 

Z = J^[t-Ml-ml^,,.. (19) 

is the on-shell wave function renormalization constant of /^^(O^) meson. 

In order to fix the unknown constant A in F{t) given by Eq. fll6p . we apply the soft-pion 
theorem 

nml) = , /"""f , (20) 
jKimc + mu) 

to the extrapolated value of the matrix element (0|c(0)s(0)|Di^') at t = m|,, from which we 
can deduce the constant 

A = „ Z^";^ . [Ml - Ml - ^{ml)] . (21) 

With all the above equipments, the DK continuum contribution to the hadronic spectral 
function can then be written as 



1 / {mD + rriKY /, {mo - rriKy 
xe{Vt - mn - ttik) 9{so - t) . (22) 



In the above calculations we have neglected the contribution of the Dgi] channel. The 
threshold of this channel is at t = 6.329 GeV^. Our formula for the contribution of the 
two-particle term is proportional to (t — Mq)"^. The lower part of the integration in t is 
more important. At the thresholds of the two channels the factor {t — Mq)"^ for the DK 
channel is about 18 times larger than that for the DgT] channel. Therefore, the effect of the 
latter is expected to be small. 
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IV. NUMERICAL RESULTS AND DISCUSSIONS 



In our numerical analysis, we use the recent result for c quark current mass rric 



rrir 



28 



1.286 GeV 
as{mz) = 0.1189 
{ss) = 0.8 X 



. Other input parameters are the following (assuming the isospin symmetry] : 



29|, m,(2 GeV) 
0.243)3 GeV^ H, {sgsa-Gs) 



96.10 MeV [30|, m„(2 GeV) = m,(2 GeV)/24.4 
0.8 GeV^x (ss) [IS^], (a.G^) = 0.06 GeV^ 
in the four quark condensation term Eq. (]A15|) a = 3, mn = "^o±+^do _ j^ggg g MeV 



rriK 



495.66 MeV 



aH, /d = 222.6 MeV [sif, fx = 159.8 MeV |3l|. 



30 

Il9|, 

3lj, 
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FIG. 2: The variation of M with when sq = 8.0 GeV^. The sohd, dashdotted, and dashed 
curves are for the case without the DK continuum contribution, g = 4.0 GeV, and g = 7.0 GeV, 
respectively. 



The renormalized coupling constant g was determined to be in the interval g = 5.1 — 
7.5 GeV in Refs. 2S, l25|. Inclusion of the contribution of DK continuum in the sum rule 
analysis of the scalar current channel will lower the g value. Since the uncertainty is large, 
we have not calculated this correction and simply allow the renormalized g to vary in the 
region g = A.O — 7.0 GeV. 

A resonance of the Z ).^ sy stem with the natural parity has has been observed experimen- 
tally at t = 8.18 GeV^ [15j. If it is an excited state of Ds{0'^), we should confine us to Sq 
smaller than and close to 8.0 GeV^. We shall first consider this case and then discuss the 
case that this resonance is not a O"*" state. The convergence of the OPE series and dominance 
of the sum by the pole and the DK continuum terms over the QCD continuum beyond sq 
constrain the Borel mass Mb in a region depending on the parameters mc{mc) and sq. Tak- 
ing rUc = 1.286 GeV, e [0.99, 2.68] GeV^ forsn = 8.0 GeV^ and M| e [0.99, 2.41] GeV^ 



for So = 7.5 GeV^. As mentioned in Refs. 23|, [33, l33|, the convergence of the perturbative 



expansion of the two-point correlation function, when written in terms of the pole quark 
mass, is rather poor, the order «<, and loop contributions being of similar size with, or 
even larger than, the leading term, while the expansion in terms of the MS running mass 
converges much faster. However, it should be noted that, even in the MS running mass 
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FIG. 3: The variation of M with when sq = 7.5 GeV^. The other captions are the same as in 
Fig. El 



scheme, the convergence of the asymptotic series in the Dg meson system is worse than the 
one found in the Bg meson system. For Ds{0^) the first order correction amounts to about 
53% and the second order to about 47% of the leading term using the values of our input 
parameters and Sq = 8.0 GeV^. The same observation has also be made in Ref. 32|. 

We first move the DK continuum contribution to the right hand side of the sum rule. 
Then we obtain the curve of M with respect to Mb by taking the derivative of the logarithm 
of both sides of the sum rule as usually done. Since this curve depends on the unknown 
parameters Mq, we have to do it self-consistently by requiring that the M value determined 
by the sum rule for the input "trial" value of Mq both lies in the middle of the stability 
window and equals roughly to Mq. For reliability of the results we also require that the ratio 
of DK contribution to the whole sum rule is not larger than about 60%. 



0.6- 



0.55- 



0.5 - 



1 . 5 



2 

fGeV^ 



2 . 5 



FIG. 4: The ratio of the DK continuum contribution as a function of with sq 
The soHd and dashed curves are for g = 4.0 GeV, and g = 7.0 GeV, respectively. 



8.0 GeV^. 



With the input mc{m^ = 1.286 GeV, we present the variation of M with M|, for sq = 
8.0 GeV^ and Sq = 7.5 GeV^ in Figs. [2] and [3|, respectively. For comparison, we also show 
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the case without DK continuum contribution with the same set of input parameters. It can 
be seen clearly from the two figures that the inclusion of the DK continuum contribution 
can lower the M value by 60 — 40 MeV. The DK continuum contributes around 45% to 60% 
of the right hand side of the final sum rule for sq = 8.0 GeV^ as shown in Fig. |H Another 
interesting point is about the vector current decay constant /o of Ds{0'^) meson. We find 
that the inclusion of DK continuum contribution lowers the decay constant /o from about 
0.185 GeV to 0.115 — 0.132 GeV for the same Sq value as can be seen from Fig. O 



0.18- 



0.16- 

> 

0.14 - 



0.12 ^ 



0.1 



1 . 5 



2.5 



FIG. 5: The vector current decay constant /o as a function of with sq = 8.0 GeV^. The other 
captions are the same as in Fig. [2j 



For mc(mj = 1.286 GeV and Sq = 7.5 - 8.0 GeV^ we found M = 2.331 ± 0.016 GeV 



being in agreement with the experimental data 2317.8 ± 0.6 MeV 3l|- For the same value 
of input parameters, we found /o = 0.128 ± 0.013 GeV, which is, however, significantly 
lower than the ones obtained in previous literature. Here we have not included the errors 
due to uncertainties in the QCD sum rule except those from the variation of the results in 
the stability window and the Sq interval, since our main interest is the central value of the 
results. The previous results already shew that the -0^(0+) mass lies in the large uncertainty 
interval of the QCD sum rule 

Now we consider the case that the new resonance found in [15J is not a O"*" state. In 
this case the sq value can only be determined by stability analysis. The results for the mass 
M found for sq = 8.5, 8.0, 7.5, 7.0 GeV^ for ^ = 7 GeV and ^ = 4 GeV are shown in Fig. 
6 and Fig. 7 respectively. The working region for sq = 8.5 GeV^ and Sq = 7.0 GeV^ are 
[0.99,2.96] GeV^ and [0.99,2.13] GeV^ respectively. The region of the Sq value can be chosen 
by requiring the least sensitivity of the results for the mass to the value of Sq. It is clear 
from these figures that this is the region between sq = 7.5 GeV^ and Sq = 8.0 GeV^ which is 
just the region chosen above for the case of a O"*" resonance at t = 8.18 GeV^ Also the best 
stabihty with respect to Mb is achieved at Sq = 7.5 GeV^. Therefore, the results obtained 
above are essentially unchanged. 

The above results show that the contribution of DK continuum, which contains the 
physics of the coupled channel effect in the formalism of QCD sum rule, is significant and 
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FIG. 6: The variation of M with when g = 7 GeV^. The soHd, dashdotted, dashed and 
dotted curves are for the case sq = 8.5 GeV^, sq = 8.0 GeV^,so = 7.5 GeV^,and sq = 7.0 GeV^ 
respectively. 
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FIG. 7: The variation of M with M| when g = A GeV^. The soHd, dashdotted, dashed and 
dotted curves are for the case sq = 8.5 GeV^, sq = 8.0 GeV^,so = 7.5 GeV^,and sq = 7.0 GeV^ 
respectively. 



is partly the reason for the unexpected low mass of 0+ cs state. Our analysis also explains 
partly why the extracted mass of the O"*" cs state from the quenched lattice QCD simulation is 
higher than the experimental value where the DK continuum contribution was not included. 
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Appendices 



APPENDIX A: THE SPECTRAL FUNCTION pqcd(s) AT THE QUARK GLUON 
LEVEL 

In this appendix, all the relevant expressions for the spectral function p{s) are given. For 



further details, we refer the readers to Ref. 23| and references therein. 



1. The perturbative spectral function 



In perturbation theory, the spectral function Pqcd(s) has an expansion in powers of the 
strong coupling constant 



Pqcb{s) = p^'\s) + p^'>{s) a{pa) + p'^Ks) aipaf + ■■■, 



.(1) 



(2)/ 



(Al) 



with a(/ia) = asipa)/'^- The leading order term p^^\s) results from a calculation of the bare 
quark-antiquark loop and is given by 



^K + m.)2.(l-^)^ 
Stt"' s 



(A2) 



and, up to order m^, the corrections in small mass can be found in Ref. 



34 



p'^Ks) 



(rric + s 2(1 — x)mcms — 2ml 



{1 + x) mcTn^ ^ (1 — 2x — a;2)m^ 
(1 — x) s 



+ 



1-x) 



(A3) 

where x = m'^/s, and the appearing quark masses correspond to the running masses in the 
MS scheme with mc{prn) and ms{prn) evaluated at the scale pm 
The order correction p^^-'(s) can be written as 



p«(.) 



167r2 



CF{mc + msf s{l- x){{l- x) 4L2(a;) + 2 Ina; ln(l - x) - (5 - 2x) 



P'r. 



xln(l-x) +(l-2x)(3-x)lnx + 3(l-3a;)ln^-^ + -(17-33x) ^, (A4) 



mt 



where L2{x) is the dilogarithmic function. The order mass corrections to the spectral 
function can be obtained by expanding the results given by [35|, l36| up to order mj, af- 
ter the higher dimensional operators have been expressed in terms of non-normal ordered 
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condensates 



Cp {thIc + nisY nicms < (1 — x) AL2{x) + 2 lnxln(l — x) — 2(4 — x) 



X ln(l -x) + 2(3 - 5x + x^) Inx + 3(2 - 3x) In ^ + 2(7 - Qs) }> , (A5) 



Cir(mc + ms) < (1 — x) 4L2(x) + 21na;ln(l — x) 



-(2 + a;)(4-x)ln(l-x) + (6 + 2x-x^)lna; + 61n^ + (8-3a;) [> , (A6) 



Pm3 ^ = -^CF{m, + m,f^^{AL2{x) + 2\nx\n{l-x) + ^— — 

OTT^ -5 I \^ ~ ^) 

^(7 + 7x-2x2), , , (Q + 7x-2x^)^ ,(2-a;2), /i^ ] 
-2^ ^ ^Inl-x +2^ ^ ^Inx + e)- In ^ 

1 — a; r iTii 



1-x) 



l-x) 



,(A7) 



-^C^(me + m,)2^<^ 2L2(x)+lnxln(l 

(13 - 24a; - 27^^ + 2x^) , , , (12 - 

^71 V2 ^"^^ ~ ^ + 

2(1 — x)^ 

(4 - 12a; + a;^ 



22x - 27x2 + 2x^ 



+ 3 



2(1 -x)3 



2(1 -x)2 

(6 - 64x + 15x2 + 11x3) 
4(1 -x)3 



Inx 



(A8) 



The three-loop, order correction p^'^\s) has been calculated by Chetyrkin and Stein- 



hauser 20|, |2l| for the case of one heavy and one massless quark. In the present analy- 
sis, we shall make use of the program Rvs.m, which contains the required expressions for 
p*^2)(s) 20, Ell- However, since the spectral function has been calculated only in the pole 
mass scheme, following Jamin and Lange 23| , in the MS scheme we still have to add to p'-^^ (s) 
the contributions resulting from rewriting the pole mass in terms of the MS mass. The two 
contributions Aip*^^) ^nd A2p^2) which arise from the leading and first order contributions, 
respectively, are given by 



Aip 



(2), 



{m^ + msYs (3-20x + 21x' 



2(1 -x)(l -3x) r. 



(2) 
m 



(A9) 



A2P 



(2)/ 



Cf (rric + rus 



sr. 



(1) 



l-x)(l-3x) 4L2(x) +21nxln(l -x) 



;i - x)(7 - 21x + Sx^) ln(l - x) + (3 - 22x + 29x2 - Sx^) Inx 



+ -(1 -x)(15-31x) 



1, 



(AlO) 
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where explicit expressions for the coefficients Vm and can be found in Appendix B. 



2. The condensate contributions 

In the following, we summarize the contributions to the spectral function Pqcd(s) coming 
from higher dimensional operators, which arise in the framework of the OPE and parameter- 
ize the appearance of non-perturbative physics. Since the spectral functions corresponding 
to the condensates contain (5-distribution contributions, we shall present directly the Borel 

^ oo 

transformed integrated quantity ull{u) = J e"'*/"pQCD(s) ds below, where u = M|, with Mb 



being the Borel mass. 

The leading order expression for the dimension-three quark condensate is well known 
with the explicit form given by 



1 _ d + ^-.j^ ^ ^l^l 



u 2mr 2u^ 



(All) 



where the expansion up to order has been included 3J] . The first order correction to the 
quark condensate can be deduced based on the fact that the mass logarithms must cancel 
once the quark condensate is expressed in terms of the non-normal ordered condensate 0, 



37 



38| with 

o r „„2 . „„2 ..2 A ^ 

ml/u 



«n^^,)(M) = ^Ci.a(m, + m,)V(gg)(r(o,^)- [i + (i-^)(in4 + ^; ^ - 

(A12) 

where r(?T,, z) is the incomplete F- function. 

The next contribution in the OPE is the dimension-four gluon condensate with the cor- 
responding expression given by 

uUf^{u) = ^(m, + m,)2(aFF)e-™'/". (A13) 

The dimension-five mixed quark gluon condensate should also be included, since it is 
enhanced by the heavy quark mass and hence still has some influence on the sum rule. 
Again the result is well known with 

nil%<,a) = - (m. + ,n.f '^^Ai^ (i _ g) ,-™?/. . (A14) 
The last condensate contribution considered in this paper is the four-quark condensate 

„n<».,(.).-.|(.-g-g).,.p, ,Ai5, 

where a is the factor representing the deviation from vacuum saturation. The contribu- 
tions of all the other higher dimensional operators are extremely small and thus have been 
neglected. 
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APPENDIX B: RELATIONSHIP BETWEEN POLE AND RUNNING MS 
QUARK MASS 



The relationship between pole and running MS quark mass is given by 23j 

m{flm) = Mpole 1 + a{fJ,a) (Aim) + a{flaf {f^a, f^m) + ■■■ 



where 



.(1) _ ^(1) 



rlnfl - 71 In 



^(2) ^ ^(2) 
m mfl 



12 + ill - Pi) ri^^o 
(1) 



(Bl) 



[B2) 



In — r + 17 (7i - Pi) In 



- ^ — 



I a 1 
7i + Pi In — 



(B3) 



The coefficients of the logarithms can be calculated from the renormalisation group |39| . 
and the constant coefficients q and r^g ^'^^ found to be 40, 41| 



,(1) 

m,0 



,(2) 
m,0 



(B4) 



Y^-yC(2)-^C(3) + 3C(2)ln2 



+ CFTn 



1111 + 1^(2) + ^^(3) --C(2)ln2 
384 2^^ ^ 8^^ ' 2^^ ^ 



- - -C(2) 
4 2^ ^ 



with 



and 



IIC^ - ATn 



71 = 2^F: 



/52 



72 



12 L 

48 



17^2 _ lOCATuf - eCpTuf 
97Ca + 9Cf - 20Tnf 



(B5) 

(B6) 
(B7) 



APPENDIX C: RELEVANT EXPRESSIONS IN THE DK CONTINUUM CON- 
TRIBUTION 

For convenience, in this appendix we collect some relevant expressions used in Sec. IIIII 
when discussing the DK continuum contribution. Firstly, we define the following loop 
integral functions Sj(t) (with t = q^) 



So(t) = 2i 



(27r)4 (p _ ^2^) [(^ _ ^)2 _ 



(CI) 



15 



2t 



kn 



(27r)4 (P - ml) [{q - kf - ml] 
1 f t + m\ 



87r2 



-f^j^ Bo{t, ml, ml) + ^ [AoimD - A,{ml)^ | , (C2) 



S2(t) 



2i 



d^k 



(27r)4 (P - [(g - k)^ - ml] 

1 r(t 



87r2 



4t 



2 \2 



Bo{t,m^,mj^) + 



t + — m|) 
it 



(C3) 



2z 



(2^)"^ (A;2 — ml) {q — k)'^ — m 



ml Bo{t, ml, ml) + Ao{ml) 



(C4) 



S4(t) 



2i 



d^k 



ekn 



(27r)4 (p - ml) [{q - k)^ - ml] 



1 { t + ml — ml 2 



87r2 



2Vt 



mj^ Bo{t,mj^,mj^) + 



K 



2Vt 



(C5) 



2i 



d^k 



2\2 



(2vr)4 (^2 



m 



K) 



{q - ky 



m 



D 



87r2 



ml BQ{t, ml, ml) + {t + ml + ml) Ao{m 



D) 



(C6) 



Here we have taken into account two intermediate states with different charges in Eqs. flCip - 
f lC6p . Apfm^) and BQ{t,mlml) is the usual one-loop scalar one- and two-point function, 
respectively [42 1 

d^k 1 



Ao[m' 



TT^ (A;2 — m^) 



m 



2 s , m 

--7i5 + ln47r) + l-ln — 
L'e /i^ 



(C7) 



Bo{t,mlml) 



d^k 



7r2 (p 



ml 



g2 - kY' 



m. 



I] 



7^; + ln47r) + In/i^ — Fo(t, m^, m: 



2> ? 



(C8) 
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where e = A — D in D-dimensional space time, /i is the introduced renormahzation scale 
regularization, and the exphcit form of the function Fo(t, m^jm^) could be 



43 



in dimensiona 
found in Ref. 

From the three linear equations for the three unknown functions fi{t) given by Eqs. f|T3|) - 
(flSj) . we can deduce the explicit expressions for the three unknown functions fi{t) 



m 

f2{t) 

with 

Y{t) 



Y{t) 
1 

W) 



Af^go [^ofK + - S0S2] {t - to) , 

2ff,go - S1S3 + S0S4 + 2 (S^ - S0S2) v^l (t - to) 



(C9) 
(CIO) 



^2/1:^0 {2E3/^ + [2E1E4 - E3(2S2 + E3) + E0E5 + 2(EiE3 - EoE4)v^]/, 



Yit) 



-E0E4 — 2E1E3E4 + E^Es + E2 (Eg — E0E5) } {t — to) ) 



(Cll) 



4 (Eo(?o' + - t) /i + 4 [(E? - E0E2) g'o + (^2 + E3 - 2Eiv^) {t - Ml) 



+ {Ml - t) UtE? - 4E4E1 + E^ - S0E5 + 4E2(E3 - Eot) - 4(EiE3 - E0E4) Vtl I 



K 



p2 
IK 



+ 



2E3E4E1 - E5E2 - SoE2 - S2 (E2 - S0E5) {Ml -t)\{t- to) 



-Acft^gl {^oIk + S0E2) {Ml - t) 



(C12) 
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With the above results, the functions A(t) and Ai(i) can be, respectively, written as 
_ i3t-ml + ml)iM',-t) ^ , . 2 ^ , {3t - + mDjM^ - t) ^ , .2 . 



, cfWoit - Mj) + (/i^g - (t - Mi){t - ml,))(t - to) 2 ^2 

+ 256/^7r^t(t - to) ^ 

, - Mq^) + (/i^^g - - - ml)){t - tp) , 

+ 256/ivr^t(t - to) ^^^""^^ 

2c/igg(t - Mq^) + + (t - Mg)(t + ml, + ml)){t - tp) , , , 

256/i7r^t(t-to) ^ ^o(m^)AoM 

(t - mf) + m|-)Ao(mf)) + (t - m|- + mf))Ao(m|-) 



6^6^ X (t-Mo^)^(m^)Ao(mJ,) 

X Bp{t,mj^,ml) - X Bp{t,ml,ml) 



8192/f7r6t 



Stt^ ' ^' 87r2(t-to) 

(t - M2)((m|, - mj^f + 2{ml + m^)t - St^) 



X Bp{t,m]j,m]i) 



32/|.7r2t 

(t + ml, - ml){cflgl{t - M^) + (/j^g - (t - Mo^)(t - ml,))(t - to)) 

256f^nH{t - tp) 

xAo{mj^)Bo{t, 171%, m\) 
{t + ml- ml){cflgl{t - M^,) + {flgl - (t - M^,){t - ml))(t - tp)) 

256/4 7r4t(t - to) 

xAp{m\)Bp{t,m\),m\) 

^ ^'^ ~ ^^^limE^t' ~ ^^^'^ ^ Mml)Ap{ml)Bp{t, ml, m\) . (C13) 

Ai(t) = iV(t)/D(t), (C14) 

N{t) = 32/^^o2 7r2(t-to) [Ao(ml,)2 + Ao(ml,)V64/^^o^'(^-^o)A(ml,)Ao(ml,) 
-32/^ tt" (t - to) [ (ml, -m\ + t) Ap{ml) - {ml -m^-t) Ap{ml) ] 
X So(t, ml,m\) - 1024/1^ gl 7rU{t- tp) Bp{t, ml, ml) , (C15) 

D{t) = 256/^ TT^ (t - to) [ (3t -ml + m^) Ap{ml) + (3t + ml- m^) Ap{ml) ] 

+(t - to) [{t-ml + ml) Ap{ml) + {t + ml- ml) Ap{ml) ] Ap{ml) Ap{ml) 

+32/1- TT^ [2cfKgl + (ml + ml + t)(t- tp) ] Ap(ml) Ap{ml) - ?>2fW { 

X [cflal + {^l to)] Ap{ml)'' + [cflgl + {ml - t) {t - tp)] Ap{ml)'' } 

+ [{ml- ml)^ - t^ ] (t - to) Ao{ml) Ao{ml) Bp{t, ml, ml) 

+32 fin' {{ml-ml + t) [c fl g'p + {ml - t) {t - tp)] Ap{ml) 

-{ml-ml-t) [cflgl + {ml -t){t-to)] Ao{ml)} Bo{t,ml,ml) 

+256/^ TT^ {4cflglt+[{ml-ml)' + 2{ml + ml)t-3t'] {t - tp)} 

xBp{t, ml, ml) + 8192/1- tt^ t (t - to) . (C16) 
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